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Abstract 

We show that the time dependent version of Sato's equation, when apphed to capacitative 
rf sheaths is no longer independent of the electric field of the space charge, and discuss the use 
of the equation for a specific sheath model. 
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1. Introduction 

A few years ago, Sato derived an expression for the current flowing in the external circuit due to the 
motion of charged particles in the gap within a discharge tube [1], for a constant voltage applied to 
the gap. Very recently, this work was generalised by Morrow & Sato [2] to include time-dependent 
applied voltages. They used an energy balance equation to obtain their results, which in both cases 
were independent of the space charge effects. 

While their methods and inferences drawn thereby are striking in their elegance and simplicity, one 
must exercise caution while applying them, for in certain situations, these can lead to results which 
are quite the reverse of the inferences drawn in their papers, as we show here. 

In a high voltage radio frequency discharge tube, there is a high concentration of positive ions at 
the electron-depleted dark space region adjoining the electrode, and a consequent build-up of field 
distortion, so that one expects a time and frequency dependence of the conduction process. The 
build-up of field distortion can therefore be modelled by a two-layer capacitor of the Maxwell- Wagner 
type [3]. 

The plasma state is characterised by equal numbers of positive and negative charges, but because of 
diffusion effects and recombinations on boundary surfaces in the discharge tube, there is charge deple- 
tion in the adjoining gas phase, resulting in the formation of a thin sheath. Electrons diffuse fastest, 
since they are lightest and have high energies, and they leave behind them a surplus of positive 
charge and a plasma potential which is positive relative to the walls. Since there is a larger number 
of charged particles in the central plasma regions of the tube, and hence better conductivity there, 
almost all of the potential drop occurs across the sheath. The clear difference in the magnitudes 
of the potential at the plasma and at the sheath, leads one to describe the plasma-sheath system 
in a discharge tube by a two-layered capacitor, with a dielectric coefficient Csh for the sheath region. 

It is shown in this work that such a representation of the plasma-sheath system leads to terms 
in Sato's equation, which depend upon the sheath potential and field. We calculate the non-zero 
contributions to the current in the external circuit from the space charge field and the sheath 
capacitance for a symmetric RF discharge for a simplified electrode geometry, in the case of the 
Godyak-Lieberman (GL) theory [4,5] for RF sheaths, using this energy balance method. 
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2. The energy balance method of Sato &c Morrow 

The continuity equations which describe the development in time t of the space charge in a gap in 
air are [6]: 

dn 

= rieaVe - HeTJVe - neUpPi - V • (rigVe - -DeVrie) 

= UeOVf, - neUp/Si - nnnpf32 - V • {ripVp - DpSJup) 

= HeVVe - nnnpP2 - V • (n„V„ - DnVUn) (1) 

where ne,np and n„ denote the electron, positive ion and negative ion densities respectively, Ve,Vp 
and v„ denote the drift velocity vectors of the electrons, positive ions & negative ions respectively, 
a is the electron ionization coefficient, r] the attachment coefficient, and /3i and 02 are the electron- 
positive ion and negative ion-positive ion recombination coefficients respectively. D^, Dp and £>„ are 
the diffusion coefficients for the electron, positive ion and negative ion respectively. 
Combining these, one obtains the equation for the net space charge density p : 

^=e^^""-''"-"-^ =-eV.r (2) 
dt dt ^ ^ 

where e is the electron charge, and the total particle flux F is given by 

r = UpVp - UnVn - UgV^ - DpVUp + DnVUn + D^VUe (3) 

One starts with these basic equations to study the electrodynamics of the charged particles in a 
discharge tube. In an RF discharge tube, since most of the potential drop in the gap occurs across 
the sheath, we treat the plasma-sheath system as a two-layer capacitor, with a dielectric constant 
Csh for the sheath and the dielectric constant of free space eo for the plasma region. 
The total current density J in the gap is 

.^er.f ,4, 

where the electric displacement D relates to the local field E through the effective complex dielectric 
constant e of the gap: 

D = eE (5) 
This can be separated into a part Dgh describing the sheath, and a plasma part Dp-. 

D = eoEp + CshEsh = Dp+ D^h (6) 

D obeys Poisson's equation : 

V • D = p = eoV • + V • {cshEsh) (7) 
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while Ep satisfies Laplace's equation: 

V-Ep = (8) 

The plasma and sheath electric fields Ep and Egh are related to their respective potential distributions 
through 

Ep = -V^pp and E^h = -^i^sh- (9) 

The energy balance equation can be used to relate the applied potential Va and the current / in the 
external circuit, to the current density and the local field in the gap through the volume integral 
Jy dv over the discharge space: 

Val= I J.Edv (10) 
Jv 

Separating out the plasma and the sheath electric fields as: 

E = Ep + E,h (11) 

and making use of (4), one obtains: 




where we have made use of the second of equations (9), performed an integration by parts, JgdS 
representing a surface integration over the closed surface of the discharge space, and then used (2). 
The second term in the right hand side of (12) can be rewritten using (7), and after again performing 
an integration by parts, this term can be written as: 

jjsh{x,t)^^^dv = j^^E,ndv + j^i,su{x,t)^^^dS (13) 

Next, we substitute (13) and (6) back into (12) to obtain: 

[er + eo^] • Epdv + -gf^Epdv - j{eT + -^)i>shdS (14) 

The first two terms within the first volume integral on the right hand side of (14) constituted 
the final form of Sato's equation derived in [2] for a time-dependent applied voltage. The boundary 
condition chosen in [2] was, that ijjah was set to zero on both electrodes, and ipp was set to zero 
on one electrode and to the applied voltage Va on the other. The authors in [2] obtained the gap 
capacitance for a simple system from the second of the two terms in the first volume integral. 
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On the other hand, by ascribing a non-trivial dielectric constant to the sheath, we have obtained 
terms additional to those obtained in [2], and these extra terms depend upon the sheath field and 
potential. 

3. Capacitative RF discharges 

Low pressure plasma chambers are widely used in the material processing industry, such as in the 
fabrication of semiconductor wafers, and in reactive ion etching. Many applications need high ion 
energies which are generated by biasing the substrate with a radio frequency (RF) current source 
[7,8]. Proper understanding of the electrodynamics involved in the sheath and near the plasma- 
sheath boundary is thus highly desirable. 

We consider a low density, low pressure plasma where a single power supply generates both the 
discharge and the RF sheath. In this high frequency regime, the sheaths are primarily capacitative 
in nature. The sheath is assumed to be collisionless. As is well known [4,5,9,10], the analysis of 
sheath models depends on the ratio of the applied RF frequency u) to the ion plasma frequency ojpi 
: when w << iVpi, the ions cross the sheath quickly and can instantaneously adjust to the applied 
field, and the properties of the sheath at different times of the RF cycle are identical to those of a dc 
sheath having a potential given by its instantaneous value (RF plus dc). However, when ui » uipi, 
the inertia of the ions prevents them from adjusting to the applied field, and they cannot respond 
to its time variation. Then the ions cross the sheath in many RF periods and they respond only to 
the dc field - their dynamics is governed by the time-averaged field in the sheath. 
In the Godyak-Lieberman (GL) theory [4,5] (which we consider now), valid in the low density, high 
frequency, high current regime [11], the ions are assumed to react only to the dc fields and not to the 
RF fields. Also it is assumed that the transit time for the ions across the sheath is large compared to 
the oscillation time. These assumptions lead the GL theory to predict monoenergetic ions impinging 
on the substrate. 

The ions are assumed to enter the sheath with a Bohm presheath velocity [12] vb ^ where 
e is the ion charge, Te denotes the electron temperature in volts and rrii is the ion mass. The 
ion sheath-plasma boundary is taken to be stationary, and it is assumed that the electrons being 
inertialess, respond to the instantaneous field. 

The GL theory holds in the regime in which the applied RF voltage is very large compared to T^, 
so that one can assume that the electron Debye length everywhere within the sheath is much 
smaller than the thickness of the ion sheath s^, implying that the electron density drops sharply 
from rie = rij at the boundary with the plasma to rig = in the sheath (at the electrode side). 
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The electron sheath penetrates into the ion sheath for a distance s{t) from the plasma- ion sheath 
boundary at a; = 0, and oscillates between a maximum thickness of Sm and a minimum thickness 
which is a few Debye lengths distant from the electrode, so that the electron sheath thickness is 
effectively Sm — s{t). 

We follow here, the analysis given in [4,5], but modified to include a finite dielectric constant for 
the sheath. 

The ion flux is conserved at the plasma-ion sheath boundary. This is expressed by: 

n^Wj = noVB (15) 

where we denote the ion density n,(a;) at the plasma-ion sheath boundary by no, and is the ion 
velocity. From energy conservation, one obtains: 

11 

= -j^iV^ - eipshix) (16) 

where ijjsh{x) is the time averaged potential within the sheath. Combining (15) and (16), one obtains 
for the ion density: 

2 - -1/2 

n,ix) = na{l - —il^shix)) (17) 

J e 

In the GL theory, a spatially uniform, sinusoidal RF current density is assumed to pass through the 
sheath: 

jRF{t) = —Jo sin cot (18) 

This current is carried by the electrons : J = —ne^, n being the electron density in the bulk 
plasma at the sheath edge. As these electrons oscillate away from the electrode, they leave behind 
a positive space charge space, pulling the positive ions there. At the electron sheath boundary, the 
displacement current given by (18) must be equated to the conduction current, for continuity : 

ds 

— eni(s)— = —Jo sin ojt (19) 
dt 

The instantaneous electric displacement Dsh{x,t) within the sheath is then given by: 

= eni{x), X > s{t) 

= 0, X < s{t) (20) 

where Dsh{x,t) = esh{x,t)Esh{x,t) . 

The time-averaged electric displacement and potential are given by: 



dx 
dx 



= e{ni{x) - ne{x)) 
= -Esh (21) 
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Tie (a;) being the time averaged electron density within the sheath. 

The electric displacement field within the sheath can be found by integrating (20): 

Dsh=e [ miOd^ (22) 

Js{t) 

This is done with the help of eqn.(19). Integrating its left hand side between the limits and s, and 
the right hand side between the limits and cot, one gets: 

e / n^{Od^ = — (1 - coswt) (23) 







Using (23) in (22), one obtains 



ni{C)d^ = e/ ni{^)d^-e ni{^)d^ , x > s{t) 
Jo Jo 

= , X < s{t) (24) 

Since the GL theory is valid for the high frequency regime, one is interested in the time-avcraged 
quantities. These can be found from s{t). Liebcrman denotes by 2<p, the phase during which x > s{t) 
: then for a; w 0, 2(t> ^ , and for x !v Sm, 2^ « 27r. Because ne{x, t) = during the part of the 
RF cycle when x > s{t), he writes: 

fieix) = (1 - '^)ni{x) (25) 

so that —4> < ujt < (j) , for x > s{t), and ut = itj) for x = s{t). 
Then, from (23) and (24), one obtains: 



1 /■+'^ 

Dsh{x) = — / Dsh{x,(jjt)d{ujt) 

^TT J_s 



= — (sin ^ — (/> COS (/>) (26) 
Prom the second of eqiis.(21) ciiid the definition of Dshi 

we have: 

Cs/t^^T^ = - — (sin cose/)) (27) 

ax WTT 

where we have assumed that the time averaging procedure allows us to factor out Csh outside the 
spatial derivative. From (19) and (17), we get 

ax So sm 

where ut was set equal to </> in (19), s to a; , and sq = j;^- 
Combining eqns.(27) and (28) we find 

esh^ = -^(1 + cos2</. - </,sin2</>)(l - ^V^.. W)'^' (29) 
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which upon integration leads to 

(1 - ^iJshix)) = 1 - A ( 3 2,/> - I cos 2</. - ^) (30) 

-le ^sh O 4 Z 

where 

L = (31) 

and wc have assumed that e^jj is independent of (j). Substituting (30) in (17), one finds the following 
expression for the ion density; 

ni = no{l (-sin20- ^cos2</.- J)} (32) 

Csh O 4 zi 

Differentiating (26) with respect to x gives us 

= A^sin0^ = ^{1 - A(3 i^2^_ ^eos2<^- ^)}"' (33) 
Since the net charge density in the sheath psh is given by 

y-Dsh = Psh , (34) 

one obtains 

Psh = (35) 

TT 

In order to calculate the sheath capacitance, one must consider the instantaneous values. 
Prom (23) and (24), the instantaneous displacement field in the sheath is 



Dsh = e / ni{C)di -{l-cosut) 

Jo ^ 



= — (coswi — cos(?i) , X > s{t) 

= , x< s{t) (36) 

Integrating both sides of (36) with respect to x, we get 

Dsh{t) = / Dsh{x,t)dx = / tshEsh{x,t)dx 

J s(t) Js(t) 

= / (cos — cos 0) sin 1 (- sin20 — — cos2(^ — — ) d(j) (37) 

^ Ju>t ^ ^sh o 4 2 J 

7riTe,„ , „ , ttL^t /15 3 5 

= — - — (3 + 4cosa;t + cos2a;i) 1 —it + -ixit + -ncosuit 

4 e^/i vl6 8 3 

1 1 5 25 \ 

+ -ujt cos 2Ljt + j^i^t COS Aijt — — sin 2a;t — — — sin Auit I (38) 
3 48 18 576 / 

where a change of variables from x to (f) has been made in (37). We now use these results of the 
GL theory modified to include a finite sheath dielectric constant €sh, to the correct form of Sato's 
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equation (14) for an RF discharge. 

The contribution to the current in the external circuit from the last integral in (14) then is : 



If it is assumed that the plasma is a good conductor, then the electric field in the plasma can be 
taken as vanishing. Then it is the last term of (39) whose contribution is relevant to the discharge 
current. 

We consider again, the GL theory. For simplicity, we consider a symmetric discharge, and circular 
electrodes of area A separated a distance d apart. Since the electric displacement and the potential 
in the sheath are assumed to vary only in the axial direction and are assumed to be uniform in the 
radial direction (implicit in the assumption that the charge density is uniform across the radius and 
varies in the axial direction only), therefore, the area can be factored out of the surface integral. If 
we make also the further assumption that esh{x,t) » esh{t) , that is, that the spatial variation of 
the sheath dielectric constant can be neglected, then one deduces, making use of (9) and (6) in (39), 
that the contribution of the term: Jg ^^f^tpshdS in (40) coming from the RF sheath to the 
gap capacitance is given by Csh , where : 



where Dgh is given in (38). Notice, that as expected, this contribution is independent of the distance 
d separating the electrodes and depends only upon the RF frequency, the electron temperature, the 
area of the electrodes, no 5 the amplitude of the current density, and the sheath dielectric constant 

In the case of RF discharges, the contribution to the total particle flux T from diffusion of the 
particles can be taken to be vanishing because the ions and electrons would be expected to react 
faster to the rapidly changing applied RF voltage, than be influenced by the diffusion gradients in 
a significant way. Thus in these cases, only the n;v; ,{l = i, e) terms in (3) would contribute to F. 
The dc contribution to the voltage across the sheath comes from the dc ion current J,. This 
contributes to the F term in (14), for the current / in the external circuit : 

2e ^^^^T 

J, = enoVB = K{ — ) ^ (41) 
where K = ^ in Lieberman's theory, while K = ^ ior Child's law. 'ipsh is given by (30): 




(39) 




(40) 




(42) 
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Performing the surface integration: 



f eniVit/jshdS 

s 



(43) 



over the surface of the discharge space, and assuming as before a siniphficd geometry of circular 
electrodes of equal area A, and a symmetric discharge, we find that this term contributes 



where Dsh is given by (38). 

In going from (29) to (30), we have made the assumption that Cgh is independent of (j). Since Cgh 
is, in fact frequency dependent: 



where r is the sheath relaxation time, it cannot be trivially factored out of the integral. One must 
use the relation: tot = (j) before performing the integration over (j). 

Discussion 

We have shown that the time-dependent version of Sato's equation is not independent of the space 
charge electric field. We consider the specific example of the Godyak-Lieberman model for a capaci- 
tative RF discharge, and show that the sheath field gives a non-negligible contribution to the gap ca- 
pacitance of the discharge tube. Space charge effects in a discharge tube can be analysed using energy 
balance methods, by modelling the system as a two-layer Maxwell- Wagner capacitor. For frequencies 
below RF frequencies also, the argument above shows that the electric field and potential of the space 
charge can give a non-zero contribution to the external current in the circuit, when the potential drop 
across the space charge is appreciable. The diffusion terms could also give a non-trivial contribution 
in these cases to the terms coming from space charge effects : Jg DVnitl)shdS , {l = i,e), D being 
the ambipolar diffusion constant. For a symmetric discharge and circular electrodes of equal area 
^, this surface term would contribute a factor: A DLno _^ |f^) ^Dgh from the ion flux, Dgh 
being given by (38). 
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